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Abstract. Emissions of halogenated very short-lived sub- 
stances (VSLS) are poorly constrained. However, their in- 
clusion in global models is required to simulate a realis- 
tic inorganic bromine (Br y ) loading in both the troposphere, 
where bromine chemistry perturbs global oxidising capac- 
ity, and in the stratosphere, where it is a major sink for 
ozone (O3). We have performed simulations using a 3-D 
chemical transport model (CTM) including three top-down 
and a single bottom-up derived emission inventory of the 
major brominated VSLS bromoform (CHBr3) and dibro- 
momethane (CH2Br2). We perform the first concerted eval- 
uation of these inventories, comparing both the magnitude 
and spatial distribution of emissions. For a quantitative eval- 
uation of each inventory, model output is compared with in- 
dependent long-term observations at National Oceanic and 
Atmospheric Administration (NOAA) ground-based stations 
and with aircraft observations made during the NSF (Na- 
tional Science Foundation) HIAPER Pole-to-Pole Observa- 
tions (HIPPO) project. For CHBr3, the mean absolute de- 
viation between model and surface observation ranges from 
0.22 (38 %) to 0.78 (115 %) parts per trillion (ppt) in the trop- 


ics, depending on emission inventory. For CH2Br2, the range 
is 0.17 (24%) to 1.25 (167%) ppt. We also use aircraft ob- 
servations made during the 2011 Stratospheric Ozone: Halo- 
gen Impacts in a Varying Atmosphere (SHIVA) campaign, 
in the tropical western Pacific. Here, the performance of the 
various inventories also varies significantly, but overall the 
CTM is able to reproduce observed CHBr3 well in the free 
troposphere using an inventory based on observed sea-to-air 
fluxes. Finally, we identify the range of uncertainty associ- 
ated with these VSLS emission inventories on stratospheric 
bromine loading due to VSLS (Bry SLS ). Our simulations 
show Br^ SLS ranges from ~ 4.0 to 8.0 ppt depending on the 
inventory. We report an optimised estimate at the lower end 
of this range (~ 4 ppt) based on combining the CHBr3 and 
CH2Br2 inventories which give best agreement with the com- 
pilation of observations in the tropics. 
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1 Introduction 

On regional to global scales, bromine (Br) chemistry plays 
an important role in atmospheric composition. In the strato- 
sphere, through coupling with analogous chlorine radicals, 
active bromine (Br x = Br + BrO) takes part in catalytic cy- 
cles (e.g. BrO-CIO) which cause large seasonal ozone (O 3 ) 
loss during polar spring (e.g. Solomon, 1999, and refer- 
ences therein). At midlatitudes, a cycle involving hydroper- 
oxyl radicals (HO 2 ) (e.g. Lary, 1996) is also significant, 
particularly during periods of elevated stratospheric aerosol 
when heterogeneous halogen activation is enhanced (Salaw- 
itch et al., 2005; Feng et al., 2007). Reduced column O 3 
increases the transmission of potentially harmful ultravio- 
let (UV) radiation to the surface, in addition to impacting 
surface temperature and climate both directly and indirectly 
(e.g. WMO, 2011, and references therein). 

In the troposphere, where understanding of halogen im- 
pacts is evolving rapidly (e.g. Saiz-Lopez and von Glasow, 
2012), Br-mediated O 3 loss is also significant (von Glasow 
et al., 2004; Yang et al., 2005), such as in the marine bound- 
ary layer (MBL) (e.g. Read et al., 2008) where biogenic 
emissions ofhalogenated species can be large (e.g. Carpenter 
and Liss, 2000; Quack and Wallace, 2003). Modelling work 
has also highlighted the importance of halogen-driven O 3 
loss in the mid-upper troposphere (Saiz-Lopez et al., 2012). 
Through reactions involving HO x (OH and HO 2 ) and NO x 
(NO and NO 2 ), bromine chemistry may indirectly perturb 
oxidising capacity and thus impact the lifetime of greenhouse 
gases (GHGs) such as methane (CH 4 ) (e.g. Lary and Toumi, 
1997). Bromine chemistry may also impact other climate- 
relevant species; e.g. bromine monoxide (BrO) is a signifi- 
cant sink for dimethyl sulfide (DMS) - a precursor for cloud 
condensation nuclei (CCN) (Breider et al., 2010). 

Sources of organic bromine include anthropogenic emis- 
sions of long-lived halons (e.g. CBrF 3 , Halon 1301) and also 
methyl bromide (CH 3 Br), whose emissions are mostly bio- 
genic (>70%) (e.g. WMO, 2011). As their production is 
regulated under the Montreal Protocol (and amendments), 
the total tropospheric bromine burden from these gases is 
now declining, from a peak observed towards the end of the 
20th century (Montzka et al., 2003). Given their long tropo- 
spheric lifetimes, these gases are a relatively minor source 
of total inorganic bromine (Br y ) below the tropopause. How- 
ever, in the stratosphere they account for ~ 75 % of the total 
Br y budget. The remainder is thought to arise from so-called 
very short-lived substances (VSLS) of predominately natural 
oceanic origin (e.g. Sturges et al., 2000; Pfeilsticker et al., 
2000). In recent years, both observational (e.g. Sioris et al., 
2006; Dorf et al., 2006, 2008; Salawitch et al., 2010; Brinck- 
mann et al., 2012) and modelling (e.g. Schofield et al., 2011; 
Hossaini et al., 2012b; Tegtmeier et al., 2012; Aschmann and 
Sinnhuber, 2013) studies have constrained their contribution 
to stratospheric Br y (Br^f SLS ) - currently estimated at 1-8 
parts per trillion (ppt) (Montzka and Reimann, 2011). 


The most abundant Br-containing VSLS are bromoform 
(CHBr 3 ) and dibromomethane (CH 2 Br 2 ) with mean MBL 
mixing ratios of ~ 1.1 and 1.5 ppt. As their nominal surface 
lifetimes are short (~26 and 120 days, assuming [OH] = 
1 x 10 6 molecules cm -3 and a global/seasonal mean photol- 
ysis rate), and their emissions exhibit significant spatial and 
temporal inhomogeneity, tropospheric gradients can be large 
(Montzka and Reimann, 2011). Localised hot-spots , where 
emissions are relatively strong, have been identified; for ex- 
ample Mace Head (Ireland) (e.g. Carpenter et al., 2005). At 
present, the total global source strength of these VSLS are 
poorly constrained and range from 430 to 1400GgBryr -1 
and 57-280 GgBryr -1 for CHBr 3 and CH 2 Br 2 (Montzka 
and Reimann, 2011). For global-scale models, a sound treat- 
ment of the magnitude and spatial distribution of VSLS emis- 
sions is required in order to simulate a reasonable Br y bud- 
get in both the troposphere and the stratosphere. As recent 
chemistry-climate model (CCM) studies suggest Br^f SLS in 
the lower stratosphere may increase in response to climate 
change (Dessens et al., 2009; Hossaini et al., 2012a), valida- 
tion of VSLS emission inventories is particularly important. 

Here, we use a three-dimensional (3-D) chemical trans- 
port model (CTM) to investigate global CHBr 3 and CH 2 Br 2 
emission inventories. We perform the first concerted evalua- 
tion of three top-down and a single bottom-up derived inven- 
tory using a combination of long-term ground-based obser- 
vations and aircraft observations. A case study for the period 
of the 2011 Stratospheric Ozone: Halogen Impacts in a Vary- 
ing Atmosphere (SHIVA) campaign, in which aircraft VSLS 
observations were obtained over the poorly sampled tropi- 
cal western Pacific, is also performed. Finally, we update 
our previous model estimate of stratospheric Br^f SLS based 
on these new emission data sets. Section 2 provides a de- 
scription of the CTM and the emission inventories. Section 3 
contains a quantitative comparison of the CTM with ground- 
based data. Section 4 contains a comparison of the CTM with 
observations from the recent HIAPER Pole-to-Pole Observa- 
tion (HIPPO) campaigns. Section 5 highlights results from 
the SHIVA western Pacific case study. Section 6 examines 
the sensitivity of Br^ SLS to emission inventories. A summary 
and conclusions are given in Sect. 7. 


2 Model and experiments 

TOMCAT is a global 3-D CTM described in Chipperfield 
(2006). The CTM runs offline and uses prescribed 6 h wind, 
temperature and humidity fields from the European Centre 
for Medium-Range Weather Forecasts (ECMWF) ERA In- 
terim reanalysis. The CTM includes a treatment of convec- 
tion, described in Stockwell and Chipperfield (1999) and fur- 
ther validated in Feng et al. (2011), based on the mass flux 
scheme of Tiedtke (1989). Vertical winds are diagnosed from 
divergence. In the boundary layer, turbulent mixing follows 
the non-local scheme of Holtslag and Boville (1993). For 
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tracer advection, the CTM uses the conservation of 2nd-order 
moments scheme of Prather (1986). The CTM was run with 
a resolution of ~2.8° longitude x ~2.8° latitude and with 
60 hybrid sigma-pressure ( a-p ) levels (surface to ~ 60 km). 

The CTM configuration here is similar to that of Hos- 
saini et al. (2012b) and includes 5 brominated very 
short-lived (VSL) source gases (SGs): CHBr3, CH2Br2, 
dibromochloromethane (CHBr 2 Cl), bromodichloromethane 
(CHBrC^) and bromochloromethane (CH2BrCl). Loss of 
these SGs occurs via oxidation with the hydroxyl radical 
(OH) or by photolysis, calculated using the recommended 
rate constants/absorption cross-section data of Sander et al. 
(2011). For simulations here, the CTM used a prescribed 
monthly mean OH field which was used in TransCom-CH4 
(Patra et al., 2011) and produced reasonable simulations of 
methyl chloroform (CH3CCI3) and CH4. 

2.1 Biogenic emissions from the ocean 

Given the significant uncertainty in global VSLS emissions, 
TOMCAT was run for the period 1 January 1997 to 31 De- 
cember 2011 with 4 previously published oceanic CHBr3 and 
CH2Br2 emission inventories. Run Sfliang used the top-down 
emission fluxes of Liang et al. (2010) (hereafter “Liang- 
2010”). Run ^Warwick used the top-down estimates described 
in Warwick et al. (2006) and updated in Pyle et al. (2011) 
(hereafter “Warwick-201 1”). Run Sbrdonez used the top-down 
estimates of Ordonez et al. (2012) (“Ordonez-2012”). Fi- 
nally, run Sziska used the bottom-up emission fluxes proposed 
by Ziska et al. (2013) (“Ziska-2013”). The global total emis- 
sions for each source gas under each scenario is given in Ta- 
ble 1. 

The top-down inventories described below rely on air- 
craft observations of CHBr3 and CH2Br2 to constrain sur- 
face fluxes using a global model. Assumptions regarding the 
latitudinal distribution and the relative importance of open 
ocean versus coastal emissions are made and vary between 
inventories. Therefore, some of the most significant uncer- 
tainty in the top-down approach is due to the lack of avail- 
able CHBr3 and CH2Br2 observations in the free troposphere 
over certain regions and due to the lack of understanding 
of emission distribution. An iterative modelling approach is 
used whereby the magnitude of emissions within a given lat- 
itude band is adjusted, in a sequential set of simulations, to 
yield the optimised agreement with observations. Model pa- 
rameters, such as coarse horizontal resolution, may add fur- 
ther uncertainty as strong local emissions, such as those from 
coastal regions, are smeared over a relatively large grid box. 
The use of observations representative of the background tro- 
pospheric CHBr3 and CH2Br2 loading may also lead to an 
under-representation of particularly strong local emissions or 
hot-spots in the top-down approach. 

The Warwick-20 1 1 scenario is a top-down estimate based 
on the original work of Warwick et al. (2006). Aircraft ob- 
servations of CHBr3 and CH2Br2, collected during the 1999 


Table 1 . Summary of 14yr CTM runs and the global total source 
strength (Gg source gas yr _1 ) of CHBr 3 and CH 2 Br 2 - 


Run 

Scenario 

Derivation 

CHBr 3 

CH 2 Br 2 

^Liang 

Liang-2010 

Top-down 

450 

62 

^Warwick 

Warwick-2011 

Top-down 

380 

113 

^Ordonez 

Ordonez-2012 

Top-down 

533 

67 

^Ziska 

Ziska-2013 

Bottom-up 

183 

64 


National Aeronautics and Space Administration (NASA) Pa- 
cific Exploratory Mission (PEM) Tropics B, were used to 
constrain surface emissions. The updated scenario used here 
is based on scenario 5 outlined in Warwick et al. (2006), how- 
ever South-east Asian CHBr3 emissions have been scaled 
down to give agreement with surface observations collected 
at Danum Valley, Borneo. This updated scenario is further 
described in Pyle et al. (201 1). 

The Liang-2010 scenario is also a model top-down es- 
timate constrained by aircraft observations. These observa- 
tions were mostly concentrated around the Pacific and North 
America between 1996 and 2008 and include the following 
campaigns: PEM-Tropics, TRACE-P, INTEX, TC4, ARC- 
TAS, STRAT, Pre-AVE and AVE (Liang et al., 2010). The 
emissions were formulated using a baseline scenario from 
Warwick et al. (2006), which was adjusted in both magni- 
tude and location, so that modelled CHBr3 and CH2Br2 gave 
good agreement with observations in the mid-troposphere, 
and the observed vertical gradient was well represented. The 
spatial distribution of emissions is assumed to be equal for 
CHBr3 and CH2Br2. 

The Ordonez-2012 scenario is the third top-down es- 
timate. It is formulated using the same aircraft observa- 
tions as Liang-2010 but also includes those obtained dur- 
ing the NASA POLARIS and SOLVE missions. This sce- 
nario is relatively sophisticated as, in the tropics (±20°), 
VSLS emissions are weighted towards the concentration of 
chlorophyll a (chi a); a potential proxy for oceanic bio- 
productivity. A monthly-varying satellite chi a climatology 
was used which allows some seasonality in the magnitude 
of the CHBr3 and CH2Br2 emission fields. This is the only 
inventory to consider such seasonality. Outside of tropical 
latitudes the sea-air flux is constant with coastal emissions 
assumed to be a factor of 2.5 larger than the open ocean. 

Finally, the Ziska-2013 scenario is a bottom-up estimate 
of emissions. Based on data of the HalOcAt database project 
(https://halocat.geomar.de/), global surface marine and at- 
mospheric concentration maps of CHBr3 and CH2Br2 (and 
CH3I) were calculated in order to derive global sea-to-air 
flux estimates. The available in situ measurements were clas- 
sified according to current knowledge about the distribution 
and possible sources of each compound, as well as the phys- 
ical and biogeochemical characteristics of ocean and atmo- 
sphere. Missing 1° x 1° grid values were extrapolated with 
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Fig. 1. Bromoform emission field (10 -13 kgm -2 s -1 ) on 1° x 1° 
grid for global (left) and western Pacific (right) regions. Emissions 
from the (a) Liang-2010, (b) Warwick-201 1, (c) Ordonez-2012 and 
(d) Ziska-2013 scenarios. 


the ordinary least square (OLS) regression technique depend- 
ing on longitude and latitude. The OLS method includes out- 
liers and thus represents the spread and variable concentra- 
tion distribution well. Based on the generated marine and 
atmospheric surface concentration maps, global climatolog- 
ical emission maps were calculated with a commonly used 
sea-to-air flux parameterisation. This applied highly tempo- 
ral (6 h) resolved wind speed, sea surface temperature, salin- 
ity and pressure data (Ziska et al., 2013). 

Global emission maps for CHBr 3 and CH 2 Br 2 are shown 
in Figs. 1 and 2, respectively. All scenarios differ signifi- 
cantly over the tropical western Pacific - an important re- 
gion for the troposphere-stratosphere transport of VSLS 
(e.g. Aschmann et al., 2009), where observations of these 
species are limited. The latitudinal-dependence of emissions 
is shown in Fig. 3. For CHBr 3 , significant variation between 
the top-down derived estimates (Warwick-20 1 1 , Liang-2010, 
Ordonez-2012) and the bottom-up estimate (Ziska-2013) is 
apparent - particularly in the tropics (±20°) and at high lati- 
tudes (> 60°) in the Northern Hemisphere (NH). For CH 2 Br 2 , 


Fig. 2. As Fig. 1 but for dibromomethane (10 14 kg m 2 s *). Note 
the change in scale. 

the total global source strength between inventories is more 
consistent (Table 1), with the exception of Warwick-20 1 1 in 
which it is ~ 1.7x larger than the others. Both Warwick- 
2011 and Ziska-2013 exhibit a significantly stronger CH 2 Br 2 
emission in the tropics relative to Liang-2010 and Ordonez- 
2012. The Ziska-2013 inventory also contains particularly 
strong emissions in the Southern Hemisphere (see Ziska 
et al., 2013). 

For minor VSLS (CHBr 2 Cl, CHBrCl 2 & CH 2 BrCl), emis- 
sions are not specified, rather their surface abundance is con- 
strained using an assumed uniform volume mixing ratio (0.3, 
0.3, 0.5 ppt) based on compiled observations in the tropical 
MBL (Montzka and Reimann, 2011). Note, Warwick et al. 
(2006) and Ordonez et al. (2012) reported emissions for these 
species but they are not available at present from the other 
inventories considered. While these minor VSLS are not the 
focus of this work, they are included in the calculation of 
Br^LS i n Sect. 6. As their nominal lifetime is relatively 
long-lived at the surface (59, 78 and 137 days) (Montzka and 
Reimann, 2011), the spatial distribution of emission is less 
important for their troposphere-stratosphere transport (rela- 
tive to CHBr 3 ). 
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Latitude (degrees) 

Fig. 3. Zonally averaged global emission source strength 
(ICr 13 kgm -2 s _1 ) for (a) CHBr 3 and (b) CH 2 Br 2 . 


3 Evaluation of emission inventories with long-term 
ground-based observations 

Previous model studies have used aircraft observations to 
validate simulated VSLS profiles in the upper troposphere 
(e.g. Liang et al., 2010; Ashfold et al., 2012; Hossaini et ah, 
2012b; Ordonez et al., 2012). Ideally, global models should 
be evaluated against observations from multiple platforms. 
For VSLS, whose emissions are poorly constrained and 
represent a significant uncertainty in global-scale models, 
a robust validation of available emission inventories with 
ground-based observations is desirable. As the troposphere- 
stratosphere transport of VSLS is highly dependant on the 
location of emissions (Aschmann et al., 2009), validation of 
both the spatial distribution and magnitude of emissions is 
needed. However, to date an evaluation of published emis- 
sion inventories has yet to be be performed. 

In this study, multi-annual observations of CHBr 3 and 
CH 2 Br 2 at 14 ground-based stations (Table 2) have been 
used to validate modelled fields and test emission esti- 
mates. The observed data are from an ongoing cooperative 
flask sampling program of the National Oceanic and At- 
mospheric Administration/Earth System Research Labora- 
tory (NOAA/ESRL). Figure 4 shows the location of obser- 
vations. Whole air samples (WAS) were collected approxi- 
mately weekly into paired steel or glass flasks and were anal- 
ysed using gas chromatography/mass spectrometry (GC-MS) 


Table 2. Summary and location of NOAA/ESRL ground-based sta- 
tions arranged from north to south. * Stations SUM, MLO and SPO 
elevated at ~ 3210, 3397 and 2810 m respectively. 


Station 

Name 

Lat 

Lon 

ALT 

Alert, NW Territories, Canada 

82.5° N 

62.3° W 

SUM* 

Summit, Greenland 

72.6° N 

38.4° W 

BRW 

Pt. Barrow, Alaska, USA 

71.3° N 

156.6° W 

MHD 

Mace Head, Ireland 

53.0° N 

10.0° W 

LEF 

Wisconsin, USA 

45.6° N 

90.2° W 

HFM 

Massachusetts, USA 

42.5° N 

72.2° W 

THD 

Trinidad Head, USA 

41.0° N 

124.0° W 

NWR 

Niwot Ridge, Colorado, USA 

40.1° N 

105.6° W 

KUM 

Cape Kumukahi, Hawaii, USA 

19.5° N 

154.8° W 

MLO* 

Mauna Loa, Hawaii, USA 

19.5° N 

155.6° W 

SMO 

Cape Matatula, American Samoa 

14.3° S 

170.6° W 

CGO 

Cape Grim, Tasmania, Australia 

40.7° S 

144. 8°E 

PSA 

Palmer Station, Antarctica 

64.6° S 

64.0° W 

SPO* 

South Pole 

90.0° S 

- 


Location of grond based and aircraft observations 

-135 -90 -45 0 45 90 135 



Fig. 4. Location of NOAA/ESRL ground-based monitoring stations. 
Note, the close proximity of stations MLO and KUM (see Table 2). 
In this work we group the stations into 5 latitude bands; high NH 
(>60°N), midlatitude NH (30-60° N), tropical (±30°), midlati- 
tude SH (30-60° S) and high SH (> 60° S). Also shown are the 
flight tracks from the NSF HIPPO aircraft campaigns (1-5) which 
took place between 2009 and 2011 (see Sect. 4). The location of the 
SHIVA aircraft campaign (see Sect. 5) that took place in the tropical 
western Pacific during Nov-Dee 2011 is also indicated. 


(Montzka et al., 2003). NOAA data from flasks collected at 
surface sites and also on the HIPPO aircraft campaign are 
presented relative to the NOAA-2003 scale for CH 2 Br 2 and 
the NOAA-2004 scale for CHBr 3 . These scales consist of 
2-4 standards prepared with gravimetric techniques at 3- 
20ppt in high-pressure (900 psi initially) 30 L, electropol- 
ished stainless steel cannisters. 

Figure 5 shows the NOAA/ESRL observed CHBr 3 mixing 
ratio at these stations (north-south). The observed data points 
are monthly mean fields that have been calculated from a 
14 yr monthly mean data set (i.e. we have taken the mean of 
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Fig. 5. Comparison of observed monthly mean CHBr 3 mixing ratio (ppt) at 14 NOAA/ESRL ground stations with output from TOMCAT 
runs ^Liang (Liang-2010 emissions), ^Warwick (Warwick-2011 emissions), ^Ordonez (Ordonez-2012 emissions) and Szi&a (Ziska-2013 emis- 
sions). The vertical bars denote ±1 standard deviation on the observed mean (see text for details). 


monthly mean fields). This approach smooths intra-monthly 
variability but can give a clear signal of seasonal variations. 
The observations spanned the period 1 January 98 to 1 Jan- 
uary 2012 at all stations except SPO, THD, and SUM, which 
are shorter records. Also shown in Fig. 5 is the correspond- 
ing modelled CHBr 3 mixing ratio from runs Su<mg, ^Warwick, 
^Ordonez and Sziska- The CTM was run for the same (14 yr) 
period following 3 yr of spin-up. Monthly mean data was out- 
put allowing a like-for-like comparison between model and 
observation. 

At NH high-latitude (> 60°) stations (ALT, SUM and 
BRW), observed CHBr 3 exhibits a pronounced seasonal cy- 
cle with elevated mixing ratios during NH winter (DJF). This 
seasonality, likely due to the enhanced photochemical sink 
of CHBr 3 during summer (JJA) months (or potentially trans- 
port), has been previously observed (at ALT) by Yokouchi 
et al. (1996). The CTM captures this seasonality, particularly 
at ALT and SUM, where the bias between model and obser- 
vation is highly dependent on the emission inventory used. 
The top-down inventories (Liang-2010, Warwick2011 and 
Ordonez-2012) on average underestimate observed CHBr 3 at 
these high-latitude NH stations (Fig. 5). The calculated mean 
bias (model minus observation) for the entire 14yr monthly- 
mean data set is —0.65, —1.61 and —0.88 ppt for these inven- 
tories, respectively. The bottom-up estimate of Ziska-2013 
overestimates with a positive mean bias of +0.54 ppt. This is 


skewed by the significant overestimation of CHBr 3 at ALT. 
It was previously shown in Fig. 3 that Ziska-2013 exhibits a 
significantly larger CHBr 3 source at high NH latitudes over 
the other inventories considered. 

At NH midlatitude (30-60° N) stations (MHD, LEF, HFM, 
THD and NWR), the agreement between model and obser- 
vation varies significantly with emission inventory. At Mace 
Head (MHD), the top-down inventories underestimate the 
large background CHBr 3 (up to ~ 8 ppt). However, the larger 
bottom-up emissions of Ziska-2013 in this region lead to a 
reasonable agreement between model and observation. Note, 
here the seasonal cycle is out of expected phase, as a CHBr 3 
minimum is observed during winter months and a maximum 
during summer. Carpenter et al. (2005) observed a similar 
seasonality and deduced that strong local emissions (during 
summer) dominate over enhanced photochemical loss to con- 
trol the local CHBr 3 abundance at MHD. 

For VSLS, transport to the stratosphere is most efficient 
in tropical regions where convection can rapidly loft bound- 
ary layer air into the mid/upper troposphere (e.g. Aschmann 
et al., 2009). At tropical (±30°) stations KUM and MLO 
there is also noticeable seasonality in observed CHBr 3 . This 
is in phase with most other NH stations and indicative of 
larger-scale processes (likely the photochemical sink) con- 
trolling the seasonality. The bias between the model and ob- 
servation is again varied and strongly dependent on emission 
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Fig. 6. As Fig. 5 but for CF^B^. 


inventory. The Ordonez-2012 emissions, which are weighted 
towards a seasonal climatology of chlorophyll a in the trop- 
ics, lead to an overestimate of CHBr 3 at each tropical sta- 
tion (KUM, MLO and SMO), and for all months. For these 
stations the mean bias is 0.12, 0.48, 0.76 and 0.07 ppt for 
runs ^Liang? ^Warwick? ^Ordonez an d 5ziska? respectively. This 
indicates that the Liang-2010 (top-down) and the Ziska-2013 
(bottom-up) derived CHBr 3 emissions perform particularly 
well at these locations in the tropical Pacific. 

In the SH, long-term observations of VSLS are particu- 
larly sparse. In the SH midlatitude (30-60° S) band, data 
from just one station is available (CGO). Here, CHBr 3 is gen- 
erally underestimated but reasonable agreement is obtained 
with the Liang-2010 and Ordonez-2012 inventories. This is 
also the case at the two high-latitude SH stations (60-90° S) 
PSA and SPO. Here, a clear seasonal cycle is apparent at 
the latter with a CHBr 3 maximum occurring during SH win- 
ter (JJA), consistent with Swanson et al. (2004) and Beyers- 
dorf et al. (2010) who note a similar seasonality. The CTM is 
able to reproduce this seasonality well, which is likely driven 
by photochemistry, and again the Liang-2010 and Ordonez- 
2012 scenarios provide the best agreement. 

For CH 2 Br 2 , a similar comparison between the observa- 
tions and the model has been performed (Fig. 6). Photolysis 
is a minor tropospheric sink for CH 2 Br 2 , which has a nomi- 
nal surface lifetime of ~ 120 days (Montzka and Reimann, 
2011), and whose dominant sink is by reaction with OH. 
As its lifetime is significantly longer than that of CHBr 3 


(~ 26 days), horizontal gradients are expected to be less pro- 
nounced. The observations show background mixing ratios 
in the range of ~ 0.5-1 .5 ppt at all stations (excluding MHD) 
with generally low variability. Seasonality is apparent at most 
sites in the NH (e.g. ALT, SUM, LEF, NWR, KUM, MLO 
etc.), and is likely due to seasonal changes to the CH 2 Br 2 + 
OH loss rate. The magnitude of relative variation is smaller 
than that for CHBr 3 due to the significantly longer lifetime 
of CH 2 Br 2 . 

The global CH 2 Br 2 source strength is relatively similar for 
3 out of the 4 inventories considered: 62-67 Ggyr -1 , among 
Liang-2010, Ordonez-2012 and Ziska-2013. However, it is 
significantly larger (113Ggyr -1 ) in the Warwick-2011 in- 
ventory. Also, the latitudinal distribution of emissions, in- 
cluding in the tropics, varies significantly between invento- 
ries (e.g. Fig. 3). At tropical stations KUM, MLO and SMO, 
CH 2 Br 2 is overestimated when using Warwick-2011 and 
Ziska-2013 emissions. At these stations, improved agree- 
ment is obtained using Ordonez-2012 and good agreement 
using Liang-2010. In the SH, between ~40 and 75° S, the 
Ziska-2013 inventory exhibits a particularly strong CH 2 Br 2 
source (see Fig. 3), not featured in the other inventories. 
Comparison of modelled CH 2 Br 2 with observations within 
this latitude range (i.e. CGO and PSA sites) show a signifi- 
cant overestimation of CH 2 Br 2 , by an approximate factor of 
2, when using the Ziska-2013 inventory. 

For a more quantitative evaluation of the modelled CHBr 3 
and CH 2 Br 2 fields with these long-term surface observations, 
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Table 3. Summary of calculated error metrics between NOAA/ESRL observed surface CHBr 3 with analogous fields from CTM runs Suang* 
^Ziska? ^Warwick and ^Ordonez- Shown is the mean bias (MB) and the mean absolute deviation (MAD) both in units ofppt. Also shown is the 
mean absolute percentage error (MAPE, see text). These fields were calculated for the entire 14yr period of available observation (1 January 
1998-1 January 2012) and for the 5 latitudinal bands shown in Fig. 4. The global values shown are a comparison for all 14 stations. 


Latitude 

Run ^Liang 

Run ^Warwick 

Run ^Ordonez 

Run S Z iska 


MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

> 60° N 

-0.65 

0.73 

33% 

-1.61 

1.61 

84% 

-0.88 

0.91 

42% 

0.54 

1.24 

79% 

30-60° N 

-0.25 

1.45 

84% 

-1.64 

1.64 

67% 

-0.35 

1.15 

59% 

-0.57 

0.77 

42% 

±30° 

0.12 

0.22 

38% 

0.48 

0.52 

77% 

0.76 

0.78 

115% 

0.07 

0.26 

36% 

30-60° S 

-0.20 

0.45 

19% 

-1.48 

1.48 

67% 

-0.54 

0.63 

26% 

-1.32 

1.32 

59% 

> 60° S 

-0.24 

0.32 

28% 

-0.75 

0.75 

61% 

-0.32 

0.40 

33% 

-0.93 

0.93 

77% 

Global 

-0.25 

0.80 

50% 

-1.04 

1.25 

72% 

-0.24 

0.87 

61% 

-0.30 

0.82 

55% 


Table 4. As Table 3 but for CE^B^. 


Latitude 

Run iSEiang 

Run ^Warwick 

Run ^Ordonez 

Run 5ziska 


MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

> 60° N 

-0.11 

0.15 

15% 

0.08 

0.14 

16% 

-0.11 

0.15 

15% 

-0.34 

0.36 

36% 

30-60° N 

-0.03 

0.28 

27% 

0.18 

0.33 

37% 

-0.00 

0.24 

24% 

0.02 

0.21 

22% 

±30° 

0.14 

0.17 

24% 

1.25 

1.25 

167% 

0.35 

0.35 

49% 

0.63 

0.63 

85% 

30-60° S 

-0.05 

0.10 

10% 

0.47 

0.49 

48% 

-0.25 

0.25 

23% 

0.93 

0.94 

92% 

> 60° S 

-0.11 

0.12 

13% 

0.43 

0.45 

55% 

-0.19 

0.19 

21% 

1.14 

1.14 

137% 

Global 

-0.02 

0.19 

20% 

0.44 

0.52 

64% 

0.00 

0.24 

27% 

0.30 

0.52 

60% 


three error metrics were calculated (Tables 3, 4); the mean 
bias (MB) (ppt), calculated using Eq.(l), the mean abso- 
lute deviation (MAD) (ppt), calculated using Eq.(2), and the 
mean absolute percentage error (MAPE), using Eq.(3), for 
the 5 latitudinal bands considered. Here, M and O denote the 
monthly modelled and observed fields for the entire 14 yr pe- 
riod of comparison, respectively. The total number of com- 
parison points ( n ) is 168. 



( 1 ) 

( 2 ) 

( 3 ) 


Based on the reported error metrics it is clear the perfor- 
mance of each inventory varies significantly by region. Fo- 
cusing on the important tropical latitude band, for CHBr 3 
the MAPE between model and observation ranges between 
36 % and 115 %. The best agreement, diagnosed by the low- 
est MAPE (36 %), is obtained from run Sziska (bottom-up 
emissions). Sprang also performs well in the tropics with a 


similar MAPE of 38 %, which is significantly lower than runs 
^Warwick and Sbrdonez- Note, small values of observed CHBr 3 
can cause large skew in the calculated MAPE (see Eq. 3). For 
CH 2 Br 2 , MAPE ranges from 24 to 166 % in the tropics. The 
best agreement is obtained from run Suang (24 %), using the 
Liang-2010 inventory which has the lowest total emissions in 
the tropics and also the lowest global total (see Fig. 3 also). 
The calculated mean bias presented in Table 4 confirms the 
significant overestimation of CH 2 Br 2 by runs ^Warwick and 
Sziska suggesting a significant overestimate of the tropical 
CH 2 Br 2 source in these emissions inventories. 

Overall, the global performance (all stations) of each sce- 
nario can be assessed based on the error metrics in Table 3. 
Globally, the best agreement between model and observation, 
for both CHBr 3 and CH 2 Br 2 , is obtained for run Suang; the 
global MAPE using the Liang-2010 inventory is 50 and 20 % 
for these species, respectively. To support this conclusion, 
further long-term observations of VSLS would be desirable, 
particularly in the tropics and in the poorly sampled South- 
ern Hemisphere. While the NOAA/ESRL observations are 
a valuable long-term record, the spatial distribution of sam- 
pling is limited in these regions. Therefore, we also consider 
recent aircraft observations of CHBr 3 and CH 2 Br 2 made dur- 
ing the HIPPO aircraft campaigns over the Pacific Basin that 
spanned global latitudes (Sect. 4). Observations of VSLS 
made in the poorly sampled tropical western Pacific during 
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the SHIVA campaign are also considered in a case study 
(Sect. 5). 

The error metrics presented in Tables 3 and 4 were com- 
puted for all months. To examine any potential systematic 
seasonal bias between the model and the observations, prob- 
ability density functions (PDF) have been computed by sea- 
son and latitude band (see supplementary material). For both 
CHBr3 and CH2Br2, no clear systematic seasonal bias is 
apparent. The skill of the model in reproducing the obser- 
vations is highly dependent on the magnitude/distribution 
of emissions, which were previously shown to vary signifi- 
cantly. The seasonality of these gases, observed at numerous 
ground-based stations, is generally well captured by asea- 
sonal emissions (3 of the 4 inventories). This suggests at 
these sites the seasonality is largely driven by photochem- 
istry; i.e. sinks that are well represented in the model. 


4 Evaluation of emission inventories with HIPPO 
aircraft data 

The HIAPER Pole-to-Pole Observations (HIPPO) project 
consisted of a series of aircraft campaigns between 2009 and 
2011 supported by the National Science Foundation (NSF). 
Five missions were conducted (January 2009, November 
2009, March/ April 2010, June 2011 and August/September 
2011). The aim of HIPPO was to make global measure- 
ments of various trace gases, including greenhouse gases 
CO2, CH4, N2O and also CO, SF6, CFCs and bromine- 
containing VSLS (Wofsy et al., 2011). Sampling spanned 
a range of latitudes, from near the North Pole to coastal 
Antarctica, on board the NSF Gulfstream V aircraft and 
from the surface to ~ 14 km over the Pacific Basin. As such, 
these comprehensive data complement the long-term obser- 
vations of VSLS discussed in Sect. 3 and allow for fur- 
ther evaluation of the model with varying emission inven- 
tories of CHBr3 and CH2Br2. The HIPPO data is particu- 
larly valuable for this analysis as it is independent, i.e. has 
not been used in constructing the emission inventories con- 
sidered. The data is archived at the following web address: 
http://www.eol.ucar.edu/projects/hippo/. 

Figures 7 and 8 compare mean observed profiles of CHBr3 
and CH2Br2, made during HIPPO 1-5, with modelled TOM- 
CAT profiles for the 5 latitude bands considered in this work. 
The observations here were collected using whole air sam- 
ples, in stainless steel and glass flasks, and analysed by 
two different laboratories by gas chromatography/mass spec- 
trometry (GC-MS); NOAA/ESRL and the University of Mi- 
ami. Mixing ratios from HIPPO are reported on the same 
calibration scale as the NOAA/ESRL ground-based station 
results. The model has here been sampled for each flight 
track to match the observations and allowing a point-by-point 
comparison throughout the profiles. To assess the skill of the 
model against the HIPPO observations, three error metrics 
were again computed; the MB (ppt), calculated using Eq.(l), 


the MAD (ppt), calculated using Eq.(2), and the MAPE, us- 
ing Eq.(3), for the 5 latitudinal bands considered. These are 
summarised in Tables 5 and 6 for CHBr3 and CH2Br2, re- 
spectively. 

In general, the model is able to reproduce the observed 
vertical distribution of CHBr3 well in all latitude bands. The 
overall skill of the model is highly dependent on the given 
emission inventory. At high latitudes in the NH (> 60° N), 
the MAPE between the model and observation ranges from 
~ 31 to 63 % across HIPPO 1-5. As was reported in Sect. 3 
using ground-based data, the best agreement, diagnosed by 
the lowest MAPE, is obtained from Skiang- We find run 
Siiang also gives the best agreement with observed CHBr3 
within the 30-60° N latitude band, where MAPE ranges be- 
tween 28 and 51 %. In the tropics (±30°), it was previously 
shown in Sect. 3 that run Sziska, with bottom-up derived 
emissions from sea-to-air fluxes, performed particularly well 
against long-term NOAA/ESRL ground-based observations. 
Based on the comparison with HIPPO observations here, run 
Sziska? which contains the lowest CHBr3 source in the trop- 
ics (e.g. Fig. 3), is again found to give the best agreement in 
the tropics. Here, the MAPE ranges from ~ 35 to 102 % with 
runs Sziska and Sordonez accounting for this lower and up- 
per limit, respectively. The significant overestimate of CHBr3 
from ^Ordonez? along with the similar reported overestimate 
found from the ground-based analysis, suggests that CHBr3 
emissions from the Ordonez-2012 inventory may be too large 
within the ±30° latitude band. 

In the SH, where the coverage provided from the ground- 
based stations is limited, HIPPO made a number of observa- 
tions. Within the 30-60° S band, the model performs reason- 
ably well with the MAPE ranging from ~ 39 to 69 %. The 
lower and upper limit is given by runs ^Warwick and Sordonez, 
respectively. At high latitudes in the SH (60-90° S), obser- 
vations are limited relative to other latitude bands, however 
some profiles are available for analysis. Note, the exception 
being for HIPPO-4, during which observations > 60° S were 
particularly sparse. Nevertheless, in this region the model 
performs reasonably well with MAPE ranging from ~ 42 to 
62 % with the best agreement from run Suang- 

To determine which CHBr3 emission inventory gives best 
agreement globally, i.e. spanning the range of latitudes cov- 
ered by HIPPO, we have also calculated a global MAPE (Ta- 
ble 5). Globally, the best agreement (lowest MAPE) between 
the model and observation is obtained from runs Suang and 
feiska? where MAPE is ~ 40 % for both. This supports the 
findings in Sect. 3, where it was also shown that the Liang- 
2010 and Ziska-2013 emission inventories give the best 
agreement with long-term NOAA/ESRL ground-based ob- 
servations of CHBr3. Note, while the global MAPE happens 
to be similar for these two runs, differences within the 5 lati- 
tudes bands are apparent. For example, in the tropics (±30°), 
as noted, the better agreement is obtained from the lower 
(Ziska-2013) bottom-up emissions (MAPE ~ 35 %). In fact, 
this is the only inventory that results in a MAPE < 50 % for 
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Table 5. Summary of calculated error metrics between CHBr 3 observed in the free troposphere during the HIPPO project (2009-2011) with 
analogous fields from CTM runs ^Liang, ^Warwick? ^Ordonez and Sziska- Shown are the MB and the MAD both in units of ppt. Also shown is 
the MAPE (see text). These fields were calculated for all observations made during HIPPO missions 1-5 for the 5 latitudinal bands shown in 
Fig. 4. A global value is also quoted for comparisons at all latitudes. 


Latitude 

Run »SLi an g 

Run 5\y arw j c k 

Run ^Ordonez 

Run .S' Zlska 


MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

> 60° N 

-0.16 

0.23 

31% 

-0.55 

0.55 

63% 

-0.03 

0.24 

37% 

-0.25 

0.26 

38% 

30-60° N 

0.04 

0.12 

28% 

-0.27 

0.29 

42% 

0.22 

0.23 

51% 

-0.25 

0.25 

42% 

±30° 

0.30 

0.31 

63% 

0.32 

0.33 

68% 

0.51 

0.51 

102% 

-0.19 

0.20 

35% 

30-60° S 

0.09 

0.13 

45% 

-0.07 

0.13 

39% 

0.19 

0.21 

69% 

-0.18 

0.19 

42% 

> 60° S 

-0.12 

0.21 

42% 

-0.37 

0.40 

60% 

0.06 

0.28 

62% 

-0.40 

0.41 

54% 

Global 

0.04 

0.20 

42% 

-0.17 

0.34 

54% 

0.21 

0.30 

65% 

-0.24 

0.25 

41% 
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Fig. 7. Comparison of observed CHBr 3 profiles (ppt) made during the NSF HIPPO project (campaigns 1-5, 2009-2011) with analogous 
modelled profiles from TOMCAT runs S^ang? S Warwick? ^Ordonez and ^ziska- All profiles shown are the average for the 5 latitudinal bands 
considered in this work and are also averaged vertically in ~ 1 km bins. The horizontal lines on the observed data denote the min-max 
variability from the mean. Note, very few observations were made during HIPPO-4 between 60 and 90° S. 
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Table 6. As Table 5 but for CF^B^. 


Latitude 

Run ^Liang 

Run ^Warwick 

Run ^Ordonez 

Run 5ziska 


MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

> 60° N 

-0.10 

0.17 

23% 

0.25 

0.29 

57% 

-0.02 

0.16 

24% 

-0.12 

0.28 

39% 

30-60° N 

-0.02 

0.15 

19% 

0.47 

0.47 

75% 

0.08 

0.17 

25% 

0.11 

0.25 

37% 

±30° 

0.12 

0.13 

16% 

1.11 

1.11 

134% 

0.24 

0.24 

29% 

0.54 

0.54 

66% 

30-60° S 

-0.01 

0.09 

13% 

0.67 

0.67 

101 % 

0.01 

0.12 

18% 

0.55 

0.55 

77% 

> 60° S 

-0.09 

0.13 

18% 

0.43 

0.43 

69% 

-0.06 

0.17 

23% 

0.55 

0.73 

96% 

Global 

-0.01 

0.13 

18% 

0.60 

0.61 

88% 

0.06 

0.17 

24% 

0.30 

0.45 

60% 
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Fig. 8. As Fig. 7 but for CH 2 Br 2 - 


CHBr 3 in the tropics, suggesting overestimated emissions in 
this region from the top-down inventories. 

The model is also able to reproduce the observed distribu- 
tion of CH 2 Br 2 well. Again, the overall skill of the model is 
highly dependent on the given emission inventory. For each 
of the 5 latitude bands considered, the best agreement be- 
tween the model and observation is obtained from run Suang- 
The calculated MAPE for this run is < 25 % within each lat- 


itude band and globally is ~ 18%. This supports the find- 
ings of the ground-based analysis presented in Sect. 3, where 
the Liang-2010 emission inventory, which has the lowest to- 
tal emissions of 62 Gg CH 2 Br 2 yr -1 (Table 1), was shown 
to perform particularly well. Note, the Ordonez-2012 inven- 
tory also performs well for CH 2 Br 2 with a global MAPE 
of ~ 24 %. This is a significantly better agreement that than 
obtained from the Warwick-2011 (87%) and Ziska-2013 
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inventories (63 %) that generally lead to overestimation of 
CH2Br2. Overall, for both CHBr3 and CH2Br2 the calculated 
biases between the model and the HIPPO aircraft data are 
consistent with, and support the findings of, the comparisons 
with the NOAA/ESRL ground-based observations. 


5 A case study in the tropical western Pacific 

The tropical western Pacific is a region of frequent and in- 
tense convection resulting in efficient transport of bound- 
ary layer air into the tropical tropopause layer (TTL) (e.g. 
Fueglistaler et al., 2009; Kruger et ah, 2009). A number of 
model studies have reported the importance of the tropical 
western Pacific for the transport of VSLS into the strato- 
sphere (e.g. Levine et al., 2007; Aschmann et al., 2009). The 
region is poorly sampled and local emissions, including those 
from farmed seaweed species (Leedham et al., 2013), are 
uncertain. Previous regional observations of VSLS include 
those made during the OP3 campaign on Borneo (Pyle et al., 
2011). Background CHBr3 was reported at ~ 1 ppt inland 
(Danum Valley) with a larger background (2-5 ppt) reported 
along the coast (Kunak). 

Figure 9 shows the modelled 2011 mean surface mixing 
ratio of CHBr3 over the tropical western Pacific. Different 
emission inventories lead to significant variation between the 
modelled CHBr3 abundance. The largest modelled CHBr3 
in this region is from Sfliang and Sbrdonez with ~3.25 and 
3.0 ppt around the northern coast of Borneo. These emis- 
sion inventories were derived with little or no observations 
in the tropical western Pacific (see Liang et al., 2010; and 
Ordonez et al., 2012). Runs ^Warwick and feiska show signif- 
icantly lower CHBr3 (~ 2 ppt) and this is likely due to the 
use of regional observations in the formulation of these in- 
ventories. Warwick-2011 was derived with regional scaling 
to give good agreement with observations made during OP3 
on Borneo, while Ziska-2013, the bottom-up estimate, in- 
cluded CHBr3 sea-air flux data measured in this region dur- 
ing the Trans-Brom cruise (Kruger and Quack, 2013). For 
surface CH2Br2 (also Fig. 9), the modelled mixing ratio is 
typically between ~ 1.0 and 1.5 ppt in the region of Bor- 
neo for all runs. The exception is run ^Warwick, where it is 
~ 1 ppt greater (i.e. ~ 2. 0-2. 5 ppt) due to the larger regional 
emissions in the Warwick-20 1 1 inventory. The remainder of 
Sect. 5 evaluates the CTM and emission inventories in this 
region using recent aircraft observations made in the free tro- 
posphere during the 2011 SHIVA campaign. 

5.1 The 2011 SHIVA campaign 

The SHIVA campaign is a European Union (EU) funded 
research project (http://shiva.iup.uni-heidelberg.de/). A pri- 
mary SHIVA objective is to investigate biogenic emissions 
of VSLS, their atmospheric transformation, transport to the 
stratosphere and ultimately their impact on O3 . A field cam- 



Fig. 9. Modelled mean surface mixing ratio (ppt) of CHBr 3 (left 
column) and CH 2 Br 2 (right column) over the tropical western 
Pacific during 2011 for CTM runs (a) ^Liang? 0*) ^Warwick? ( c ) 
^Ordonez an d ( A ) ^Ziska- 


paign was conducted during November-December 2011 in 
the tropical western Pacific region based on Malaysian Bor- 
neo. An overview of the campaign is given in Pfeilsticker and 
the SHIVA consortium (2013). 

5.1.1 Aircraft observations 

Aircraft observations of VSLS in the tropical western Pa- 
cific region are extremely limited. Within the framework 
of SHIVA, aircraft observations of brominated VSLS were 
made during 14 flights on board the Deutschen Zentrums 
fur Luft- und Raumfahrt (DLR) Falcon aircraft around Bor- 
neo. The flight tracks and location of sampling is shown 
in Fig. 10. Here we consider observations of major VSLS 
CHBr3 and CH2Br2 made by the University of Frankfurt 
(UOF) and the University of East Anglia (UEA). These data 
are used to further evaluate the performance of the model, 
and top-down/bottom-up emission inventories, in the free 
troposphere within this poorly sampled region. 

Observations made by the UOF group used the Gas chro- 
matograph for Observation of Stratospheric Tracers-Mass 
Spectrometer (GhOST-MS) instrument - a fully automated 
GC/MS system for airborne (in situ) observations of halo- 
genated hydrocarbons. Observed mixing ratios for CHBr3 
and CH2Br2 from the GhOST-MS are reported on the 
NOAA-2003 calibration scale (see Sects. 3 and 4). The de- 
termined accuracy of the working standard gas is estimated 
at 16.5 and 9.0% for these species, respectively. The preci- 
sion of the instrument varies between flights but is typically 
< 4 % for both species. For further details of the SHIVA air- 
craft observations see Sala et al. (2013). 

Observations by UEA used the Falcon’s whole air sam- 
pler (WASP) that consisted of 30 glass flasks (approximately 
700 mL internal volume) which were filled to a pressure of 
2.5 Bar using a diaphragm pump. The samples were analysed 
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Fig. 10. Flight tracks of the DLR Falcon aircraft during November and December 2011 as part of the 2011 SHIVA campaign. 


for halocarbons within 48 h of collection using a GC/MS 
(Agilent 6973) operating in negative ion, chemical ionisation 
mode (Worton et al., 2008). Because of a limitation of the 
sampling pump, WASP samples were only collected at alti- 
tudes below ~ 3 km. WASP data for CHBr 3 and CH 2 Br 2 are 
also reported on the most recent NOAA scales. Typical ana- 
lytical precision (750 mL sample) was <4% for both com- 
pounds, with a calibration uncertainty of 7.1 and 6.5% for 
CHBr 3 and CH 2 Br 2 , respectively. The two bromocarbon data 
sets will be examined in further detail in Sala et al. (2013). 

Figure 11 shows the modelled mixing ratio of CHBr 3 
sampled along the flight track of the Falcon aircraft dur- 
ing SHIVA. Also shown is the observed CHBr 3 mixing ra- 
tio from the GhOST-MS and WASP instruments. The ob- 
servations show that during most flights, CHBr 3 rarely ex- 
ceed 1.0-1.5ppt. A notable exception is flight 4, during 
which CHBr 3 was elevated (>2ppt) near the surface. Large 
quantities of seaweed were visible from the aircraft during 
this flight, suggesting a large and localised emission source. 
Note, within the framework of SHIVA, emissions of halo- 
carbons from both naturally occurring, and farmed tropical 
macroalgae, has been investigated (Leedham et al., 2013). 

The agreement between modelled and observed CHBr 3 is 
highly dependent on the emission inventory. As before, we 
have calculated the MB (ppt), the MAD (ppt) and the MAPE 
(%) between the model and observation for all flights con- 
sidered. For CHBr 3 , a summary of these metrics is given in 
Table 7. In general, the top-down inventories (Liang-2010, 
Warwick-2011, and Ordonez-2012) overestimate the obser- 
vations. This is particularly the case for runs Suang and 
^Ordonez where CHBr 3 is overestimated, from the surface up 
to ~ 12 km, during numerous flights (e.g. flights 2a, 7, 10b). 
The MB between model and observation for these flights is 
1.67, 1.32 and 0.96 ppt for Suang and 1.61, 1.19and 0.99 ppt 
for ^Ordonez- Whilst also overestimating, an improved agree- 
ment is obtained from run ^Warwick in this region. For exam- 
ple, for the same flights the MB is smaller (i.e. MB < lppt) 
at 0.82, 0.78, and 0.47 ppt. Overall, the best agreement is ob- 
tained from Sziska (bottom-up emissions), which for some 


flights exhibits a small negative bias. For the above flights, 
the MB from this run is 0.31, 0.06 and 0.07 ppt, respectively. 

Across all the flights considered, the MAPE between the 
model and observed CHBr 3 is 117, 68, 125 and 37% for 
runs ^Liangj ^Warwick? ^Ordonez and 5ziska? respectively — high- 
lighting the significant variation in the performance of the 
inventories in this region. The bottom-up CHBr 3 emissions 
proposed by Ziska et al. (2013) perform particularly well as 
this is the only inventory that gives rise to a MAPE <50% 
in this region. This inventory was also shown to perform 
well against the NOAA/ESRL ground-based observations 
(Sect. 3) and HIPPO aircraft observations (Sect. 4) in the 
tropical Pacific Basin. The Ziska-2013 inventory is con- 
strained by local sea-to-air fluxes obtained in the tropical 
western Pacific during ship cruises; e.g. Trans-Brom (Kruger 
and Quack, 2013). This is the likely explanation as to why 
the MAPE is significantly lower for this inventory, over 
Liang-2010 and Ordonez-2012, that are based on limited or 
no regional (aircraft) observations. The same is true of the 
Warwick-201 1 inventory, which also performs relatively well 
in this region, and is constrained by local (ground-based) 
observations. This further highlights the need for more lo- 
cal observations of VSLS, particularly in poorly sampled re- 
gions, in order to improve VSLS emission inventories at the 
regional scale. 

Figure 12 shows the modelled versus observed CH 2 Br 2 
during SHIVA flights. The observations show CH 2 Br 2 typ- 
ically in the range of 0.5-1. 5 ppt during most flights and 
with a relatively small vertical gradient. The performance of 
each emission inventory is assessed using the error metrics 
summarised in Table 8. Across all flights, the MAPE is rela- 
tively low (compared with that for CHBr 3 ) at 25, 119, 34 and 
56 % for runs ^Liang? ^Warwick? ^Ordonez an d Sz iska? respec- 
tively. Consistent with the NOAA/ESRL ground-based anal- 
ysis (Sect. 3) and also the HIPPO aircraft analysis (Sect. 4), 
the best agreement (diagnosed by lowest MAPE) between 
modelled and observed CH 2 Br 2 , is obtained by Suang- For 
run ^Warwick, which was previously shown to overestimate 
surface CH 2 Br 2 at NOAA/ESRL stations in the Pacific Basin 
(e.g. Fig. 6), we again find an overestimate against SHIVA 
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Table 7. Summary of calculated error metrics between CHBr 3 observed in the free troposphere during 14 flights of the SHIVA aircraft 
campaign (November-December 2011) with analogous fields from CTM runs Suang* ^Warwick? ^Ordonez an d Sziska- Shown is the MB and 
the MAD both in units of ppt. Also shown is the MAPE (see text). These fields were calculated for all observations from both instruments 
deployed during SHIVA (i.e. GhOST-MS and WASP, see text). A mean value for all 14 flights is also reported. 


Flight 

Run »SLi an g 

Run ^Warwick 

Run ^Ordonez 

Run Sziska 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

2a 

1.67 

1.67 

192% 

0.82 

0.82 

99% 

1.61 

1.61 

187% 

0.31 

0.39 

46% 

2b 

0.58 

0.58 

101% 

0.43 

0.43 

77% 

0.70 

0.70 

122% 

-0.08 

0.16 

24% 

3 

0.49 

0.61 

40% 

-0.20 

0.34 

22% 

0.55 

0.62 

44% 

0.34 

0.58 

35% 

4 

0.70 

0.70 

82% 

0.25 

0.31 

42% 

0.84 

0.85 

101% 

0.13 

0.27 

32% 

5 

0.43 

0.53 

67% 

0.04 

0.41 

51% 

0.55 

0.63 

80% 

-0.76 

0.79 

48% 

6a 

0.42 

0.57 

122% 

0.18 

0.42 

84% 

0.48 

0.59 

127% 

-0.28 

0.35 

39% 

6b 

0.62 

0.75 

113% 

0.26 

0.43 

65% 

0.73 

0.87 

127% 

-0.07 

0.27 

28% 

7 

1.32 

1.32 

308 % 

0.78 

0.78 

186% 

1.19 

1.19 

277% 

0.06 

0.20 

39% 

8b 

0.81 

0.81 

120% 

0.33 

0.35 

62% 

0.86 

0.86 

128% 

-0.08 

0.20 

28% 

9 

0.56 

0.59 

91% 

0.36 

0.43 

68% 

0.70 

0.71 

109% 

-0.23 

0.28 

34% 

10a 

1.01 

1.01 

150% 

0.40 

0.41 

74% 

1.03 

1.03 

155% 

0.18 

0.29 

40% 

10b 

0.96 

0.96 

172% 

0.47 

0.49 

106% 

0.99 

0.99 

178% 

0.07 

0.25 

40% 

11a 

0.45 

0.57 

77% 

-0.02 

0.43 

47% 

0.65 

0.71 

92% 

-0.50 

0.54 

37% 

lib 

0.49 

0.59 

78% 

-0.01 

0.38 

43% 

0.70 

0.74 

95% 

-0.42 

0.46 

34% 

All 

0.76 

0.81 

117% 

0.28 

0.46 

68% 

0.84 

0.87 

125% 

-0.12 

0.39 

37% 
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Fig. 11. Comparison between modelled and observed CHBr 3 mixing ratio (ppt) along the flight tracks of the DLR Falcon aircraft during the 
2011 SHIVA campaign. Model output is from CTM runs ^Liang » ^Warwick? ^Ordonez and ^ziska and observed data from the GhOST in situ 
GC/MS system and the WASP whole air sampler (Sala et al., 2013). The dashed pink line denotes the altitude of the aircraft. 
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Table 8. As Table 7 but for 


Flight 

Run ^Liang 

Run 5\\r arw i c k 

Run ^Ordonez 

Run Sziska 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

MB 

MAD 

MAPE 

2a 

0.38 

0.42 

45% 

1.23 

1.23 

134% 

0.45 

0.48 

52% 

0.64 

0.65 

70% 

2b 

0.16 

0.16 

20% 

1.14 

1.14 

133% 

0.27 

0.27 

32% 

0.49 

0.49 

58% 

3 

0.12 

0.13 

13% 

1.08 

1.08 

102% 

0.23 

0.23 

23% 

0.66 

0.66 

62% 

4 

0.22 

0.22 

26% 

1.08 

1.08 

125 % 

0.34 

0.34 

40% 

0.57 

0.57 

65% 

5 

0.04 

0.12 

12% 

0.90 

0.90 

89% 

0.15 

0.17 

18% 

0.22 

0.24 

25% 

6a 

0.13 

0.20 

25% 

1.00 

1.00 

124% 

0.21 

0.25 

31% 

0.37 

0.38 

49% 

6b 

0.17 

0.17 

21% 

0.97 

0.97 

115% 

0.27 

0.27 

32% 

0.41 

0.41 

49% 

7 

0.46 

0.46 

62% 

1.49 

1.49 

200% 

0.50 

0.50 

68% 

0.75 

0.75 

101% 

8b 

0.20 

0.20 

24% 

1.07 

1.07 

125 % 

0.30 

0.30 

35% 

0.52 

0.52 

61% 

9 

0.12 

0.13 

15% 

1.13 

1.13 

125 % 

0.25 

0.25 

29% 

0.59 

0.59 

65% 

10a 

0.26 

0.26 

29% 

1.10 

1.10 

126% 

0.35 

0.35 

40% 

0.57 

0.57 

65% 

10b 

0.25 

0.26 

32% 

1.18 

1.18 

139% 

0.35 

0.35 

42% 

0.65 

0.65 

77% 

11a 

0.05 

0.18 

18% 

0.88 

0.88 

94% 

0.17 

0.21 

23% 

0.27 

0.30 

34% 

lib 

0.06 

0.13 

14% 

0.89 

0.89 

92% 

0.19 

0.19 

21% 

0.30 

0.31 

33% 

All 

0.18 

0.22 

25% 

1.07 

1.07 

119% 

0.28 

0.29 

34% 

0.49 

0.49 

56% 
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Fig. 12. As Fig. 11 but for CF^B^. 


observations (approximate factor of 2). Therefore, it seems 
highly likely that the CH 2 Br 2 emission strength is signifi- 
cantly overestimated in the tropics by the Warwick-20 1 1 in- 
ventory. 


6 Sensitivity of stratospheric bromine loading to 
emission inventory 

In our previous modelling work, emissions of major VSLS 
CHBr 3 and CH 2 Br 2 were not specified in the TOMCAT 
CTM (Hossaini et al., 2010, 2012b). Rather, a uniform sur- 
face mixing ratio (~ 1 .2 ppt) was imposed in the tropics 
(±20°) based on compiled aircraft observations. Using this 
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approach, Hossaini et al. (2012b) quantified stratospheric 
Biy SLS as ~ 5 ppt; i.e. within the compiled range of 1-8 ppt 
outlined in WMO (2011), and in general agreement with 
balloon-borne estimates (Dorf et al., 2006, 2008). The CTM 
performed reasonably well against aircraft observations in 
the TTL. However, this approach meant regional hot-spots , 
where emissions may be large and background concentra- 
tions elevated, were not captured. Any dependence of strato- 
spheric Br^ SLS on the spatial distribution of surface emis- 
sions was also not modelled. Here, using the CTM runs pre- 
sented in this paper (i.e. multiple emission inventories for 
CHBr 3 and CH 2 Br 2 ), we revise our estimate of Br^ SLS based 
on these spatially varying, and seasonally varying in the case 
of Ordonez-2012, emission inventories. 

It is thought that VSLS contribute to the stratospheric 
bromine budget via both source gas injection (SGI) and 
also product gas injection (PGI). The SGI pathway is quan- 
tified by summing the total organic bromine from VSLS 
reaching the lower stratosphere. For PGI, which refers to 
the troposphere-stratosphere transport of inorganic product 
gases (e.g. BrO, HBr), the tropospheric partitioning of Br y 
among soluble and non-soluble species needs consideration. 
As this involves complex heterogeneous and multi-phase 
processes (e.g. Aschmann and Sinnhuber, 2013), which are 
crudely treated in global models, Br y speciation and recy- 
cling represents a significant uncertainty in the quantifica- 
tion of PGI with models. The approach used here is identical 
to that described in Hossaini et al. (2012b). Once Br y is re- 
leased from source gases it is partitioned between soluble and 
non- soluble form according to a mean altitude-dependent 
HBr : Br y ratio. This was taken from a previous CTM in- 
tegration in which detailed partitioning of tropospheric Br y 
was considered. 

Figure 13 shows the modelled tropical mean profile of 
Biy SLS in the stratosphere at the end of the 14yr simulation. 
We find Br^ SLS ranges from ~ 5 to 8 ppt (above ~ 30 km) de- 
pending on the choice of emission inventory. Runs Sziska and 
^Warwick account for the lower limit and upper limit, re- 
spectively. However, as ^Warwick overestimated both CHBr 3 
and CH 2 Br 2 significantly in the tropics, it seems likely that 
the upper limit of ~ 8 ppt reported here is also an overes- 
timate. We have therefore now identified a range of uncer- 
tainty with regard to emissions of major VSLS CHBr 3 and 
CH 2 Br 2 on stratospheric Br^f SLS loading. Note, here Br^f SLS 
also includes the contribution from minor VSLS CHB^Cl, 
CHBrCl 2 and CH 2 BrCl. Their total contribution to Br^ SLS is 
~ 1 ppt and is consistent between each model run. 

The modelled stratospheric Bry SLS ranges from to 
8 ppt when both CHBr 3 and CH 2 Br 2 are taken from the same 
inventory. However, in the tropics, where the troposphere- 
stratosphere transport of VSLS is most rapid, it was shown 
using ground-based (Sect. 3) and aircraft (Sect. 4) observa- 
tions that a single inventory does not provide the simultane- 
ous best agreement for both VSLS in this region. For CHBr 3 , 


Modelled tropical mean Br y (VSLS) 



Fig. 13. Modelled 2011 tropical (± 30 °) mean profile of total in- 
organic bromine (ppt) from CHBr3, CH2Br2, CHB^Cl, CTLBrCl 
and CHBrCl2 (Br y SLS ) in the stratosphere. Profiles are shown for 
CTM runs S Liang , % a rwick> ^Ordonez, and S Z iska- An optimised es- 
timate, calculated by combining CHBr3 from Sziska and CH2Br2 
from ^Liang? * s al so shown. 

the best agreement was obtained from run Sziska and simi- 
larly, for CH 2 Br 2 , run Suang gave the best agreement. There- 
fore, we also report an optimised estimate of stratospheric 
Br^ SLS based on a combination of these two fields; ~4ppt 
(also shown in Fig. 13). From the 4 inventories considered, 
the CHBr 3 and CH 2 Br 2 source strengths are the lowest in 
Ziska-2013 and Liang-2010, respectively. Therefore, the op- 
timised estimate reported here is lower than the range ob- 
tained when considering emissions of both species from the 
same inventory. 

Our optimised Br^ SLS estimate of ~4ppt is lower than 
that reported in our previous work (~ 5 ppt) (Hossaini et al., 
2012b), which did not use spatially varying emission fluxes. 
The use of a fixed mixing ratio as a surface boundary con- 
dition for CHBr 3 and CH 2 Br 2 in Hossaini et al. (2012b) 
may have overestimated their abundance in the boundary 
layer. However, our modelled optimised estimate is in good 
agreement with Br^f SLS derived from observations of strato- 
spheric BrO (the so-called inorganic method). For example, 
using differential optical absorption spectroscopy (DOAS) 
to obtain BrO profiles, combined with photochemical mod- 
elling, Dorf et al. (2006) reported a Br^ SLS contribution 
of 4.1 (±2.5) ppt. However, given the recent findings of Kr- 
eycy et al. (2013) on the ratio of 7(Br0N02)/k(Br0+N02), 
this estimate may need to be revised downward. Overall, 
our model calculations are consistent with the broad Br^ SLS 
range of 1-8 ppt reported by WMO (2011). 
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7 Summary and conclusions 

Global models require a realistic treatment of biogenic 
bromine emissions in order to simulate a reasonable Br y 
budget in both the troposphere and the stratosphere. At 
present, oceanic emissions of brominated VSLS are poorly 
constrained and represent a significant uncertainty in global 
models (WMO, 2007, 2011). Given suggestions that strato- 
spheric VSLS loading may increase in response to climate 
change (Dessens et al., 2009; Hossaini et al., 2012a), con- 
straining both the magnitude and spatial distribution of con- 
temporary emissions is important. In this study we have used 
a global model to perform the first concerted evaluation of 
previously published global CHBr 3 and CH 2 Br 2 emission in- 
ventories. We have evaluated three top-down and a bottom- 
up derived inventory by comparing the simulated abundance 
of these VSLS with independent observations - i.e. the ob- 
served data was not included in the formulation of the emis- 
sion inventories. The observed data have included long-term 
observations at various NOAA/ESRL ground-based stations, 
aircraft observations made during the NSF HIPPO cam- 
paigns (1-5) and also novel aircraft observations made dur- 
ing the 2011 SHIVA campaign over the poorly sampled tropi- 
cal western Pacific. We have also updated our previous model 
estimate of Br^ SLS based on these available emission scenar- 
ios. 

Our comparisons reveal the TOMCAT CTM is able to re- 
produce a variety of global CHBr 3 and CH 2 Br 2 observations. 
The agreement between the model and the observation is 
highly dependent on the choice of emission inventory, which 
differ significantly in terms of magnitude and spatial distribu- 
tion. All the inventories considered give good agreement in 
some locations. However, to accurately diagnose the source 
gas injection of VSLS into the stratosphere, simulating their 
abundance in the tropics, where transport to the stratosphere 
is rapid, is most important. Comparison of the model with 
observations at NOAA/ESRL surface sites and also with air- 
craft observations obtained during HIPPO, shows a consis- 
tent pattern on the performance of individual emission in- 
ventories. Based on these comparisons, along with the results 
from the SHIVA case study, our main findings are the follow- 
ing. 

- Current global emission inventories of CHBr 3 and 
CH 2 Br 2 , which are used in global models, vary sig- 
nificantly. Evaluating these inventories is challenging 
due to the limited spatial coverage of long-term obser- 
vations, particularly in the tropics and in the Southern 
Hemisphere. Averaged globally, the best agreement 
between modelled CHBr 3 and CH 2 Br 2 with long- 
term surface observations made by NOAA/ESRL is 
obtained using the top-down emissions proposed by 
Liang et al. (2010). Globally, the mean absolute per- 
centage error between the model and NOAA/ESRL 
observations for this inventory is ~ 50 and ~ 20 % for 


CHBr 3 and CH 2 Br 2 over a 14yr period, respectively. 
Comparison of the model with aircraft observations 
made during the HIPPO project, which spanned global 
latitudes over the Pacific Basin, also support these find- 
ings. Globally, the mean absolute percentage error be- 
tween the model and HIPPO observations is similar 
at 42 and 18%, for CHBr 3 and CH 2 Br 2 respectively, 
when using the Liang et al. (2010) emissions. Glob- 
ally, we also find the CH 2 Br 2 emissions of Ordonez 
et al. (2012) perform particularly well with a mean ab- 
solute percentage error of less than ~30% between 
model and observations. 

- For CHBr 3 , within the tropics only, the best agreement 
between the model and observations is obtained us- 
ing the bottom-up emission fluxes proposed by Ziska 
et al. (2013). Using this inventory, the mean abso- 
lute percentage error between the model and long-term 
NOAA/ESRL surface observations is ~ 36 %. Against 
the HIPPO observations it is ~ 35 %, with the other in- 
ventories considered giving a significantly larger bias 
(> 60 %). For CH 2 Br 2 , in the tropics the model is able 
to reproduce observations well using the Liang et al. 
(2010) emissions. The mean absolute percentage error 
is 24 and 16 % when the model is compared with the 
NOAA/ESRL and HIPPO observations, respectively. 
Therefore, our results show the bias between model 
and observations is lowest when using the inventories 
that have the lowest tropical (and also total global) 
CHBr 3 and CH 2 Br 2 emission strength (i.e. Ziska and 
Liang, respectively). 

- In the tropical western Pacific, where rapid 
troposphere-stratosphere transport coincides with 
relatively large VSLS emissions, the model is able 
to reproduce novel aircraft observations of CHBr 3 
and CH 2 Br 2 made during the 2011 SHIVA campaign. 
The skill of the model is highly dependant on the 
choice of emission inventory. Good agreement is 
obtained for CHBr 3 using the bottom-up emissions 
of Ziska et al. (2013). These emissions, along with 
the Warwick-2011 inventory, were formulated using 
local observations around Borneo and as such perform 
significantly better than inventories based on limited 
aircraft observations alone. Therefore, this further 
highlights the need for more observations of VSLS, 
particularly at a finer spatial resolution, in order to 
improve current regional emission estimates. For 
CH 2 Br 2 , which is a longer-lived source gas, the Liang 
et al. (2010) emissions were again found to give the 
best agreement with the observations. 

- The modelled contribution of VSLS to stratospheric 
inorganic bromine varies significantly depending on 
the CHBr 3 and CH 2 Br 2 inventory used. We find 
Bry SLS ranges from ~5.0 to 8.0 ppt when using 
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CHBr 3 and CH 2 Br 2 from the same inventory. How- 
ever, we find no single inventory provides the simulta- 
neous best agreement with observations in the tropics. 
Therefore, we also report an optimised estimate, calcu- 
lated by combining the inventories which perform the 
best in this region. A combination of CHBr 3 emissions 
from Ziska et al. (2013) and CH 2 Br 2 emissions from 
Liang et al. (2010) lead to our (lower) optimised esti- 
mate of — 4 ppt. These inventories were found to con- 
sistently perform the best in the tropics using three in- 
dependent sets of observations (i.e. NOAA/ESRL sur- 
face, HIPPO aircraft and SHIVA aircraft data). Both 
the modelled range and optimised estimate are within 
the compiled 1-8 ppt range reported by WMO (201 1). 
Therefore, in this study we have now identified the 
range of uncertainty associated with emissions of ma- 
jor VSLS CHBr 3 and CH 2 Br 2 on stratospheric Br^ SLS 
loading. Although, model estimates of the product gas 
injection contribution to Br^ SLS remain a significant 
uncertainty (e.g, Salawitch et al., 2010; Aschmann and 
Sinnhuber, 2013). 

Our study has shown that in recent years understanding 
of oceanic VSLS emissions has improved significantly and 
that current inventories used in global models are reason- 
able. Based on the results of this work, it would be useful to 
revise current inventories and/or combine them to improve 
the treatment of CHBr 3 and CH 2 Br 2 emissions in global 
models. Furthermore, it would be useful for the analysis 
performed in this study to be repeated by other modelling 
groups, in order to determine the extent to which our re- 
sults are model dependent. For example, to assess the role 
of differences in model transport, such as mixing in the 
boundary layer and convection, which are parameterised and 
likely to vary between models. A related exercise examining 
inter-model variability on the tropospheric distribution and 
the troposphere-stratosphere transport of VSLS is planned, 
within the framework of the ongoing Atmospheric Tracer 
Transport Model Inter comparison (TransCom) project (e.g. 
Patra et al., 2011). This project, TransCom- VSLS, will ex- 
amine the performance of a number of global models against 
existing VSLS observations and also assess the variability 
between data sets, including the impact of temporal sampling 
and systematic biases on the agreement between models and 
observations. Finally, future work will examine emissions of 
relatively minor VSLS (e.g. CHB^Cl, CH 2 BrCl) along with 
a more detailed examination of emission seasonality. 
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